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NUCLEAR  RESONANCE  PROFILING  DATA  ANALYSIS 


I.  INTRODUCTION 

The  purpose  of  this  report  is  to  describe  some  of  the  methods  devel¬ 
oped  at  NHL  during  the  past  decade  for  nuclear  resonance  profiling  (NRP) 
with  the  principal  objective  of  providing  aid  to  those  at  NRL  who  wish  to 
use  these  methods  for  the  acquisition  and  analysis  of  pertinent  data. 

In  this  report,  NRP  will  mean  the  acquisition  and  analysis  of  nuclear 
reaction  data  obtained  with  incident  ion  energies  in  the  neighborhood  of  a 
narrow  nuclear  resonance  for  the  purpose  of  measuring  the  absolute  concen¬ 
tration  of  some  element  as  a  function  of  depth  in  the  near-surface  region  of 
a  solid.  The  term  near-surface  is  meant  to  designate  a  region  a  micrometer 
or  less  below  a  well-defined  surface  of  a  solid  with  cross-sectional  area 
the  size  of  the  incident  ion  beam.  Concentration  as  a  function  of  depth 
will  be  referred  to  as  a  concentration  function  or  a  profile. 

The  methods  described  herein  are  rooted  in  the  development  at  NRL  of 
a  versatile  Van  de  Graaff  accelerator  with  a  set  of  sophisticated  ancillary 
apparatus,  research  in  nuclear  physics  (especially  in  connection  with  sharp 
(p t'f)  resonances),  and  the  need  for  precise  "impurity  profiles"  in  the  near¬ 
surface  region  of  materials  useful  in  the  development  and  fabrication  of 
solid-state  electronic  devices.  Emphasis  in  this  development  was  given  to 
situations  where  existing  methods  of  measurement  were  inadequate. 

’■Then  the  use  of  NRP  methods  is  contemplated,  a  large  collection  of 
methods  based  upon  the  bombardment  of  a  specimen  with  particles  or  electro¬ 
magnetic  radiation  and  the  consequent  emission  of  particles  and/or  radiation 
should  be  considered.  An  impressive  array  of  such  methods  has  been  devel¬ 
oped  within  the  past  decade  and  most  of  these  are  more  convenient,  more  rapid, 
and  less  expensive  than  NRP.  However,  there  are  problems  in  material  analy¬ 
sis  where  NRP  can  yield  results  not  obtainable  by  other  means.  For  example, 
much  of  the  impetus  for  NRP  at  NRL  came  from  the  need  for  sodium  impurity 
profiles  in  Si02  &bout  1972.  It  was  found  that  sodium  distributions  changed 
when  measurements  were  attempted  by  means  of  secondary-ion  mass  spectroscopy 
(SIMS)  but  not  during  NRP  measurements  (1). 

Detection  limit  is  an  important  parameter  in  materials  analysis  and  is 
here  taken  to  mean  the  minimum  concentration  (of  the  element  of  interest  in 
the  test  specimen)  for  which  useful  information  can  be  obtained.  A  particular 
detection  limit  must  be  associated  with  a  particular  set  of  measurement 
parameters.  Vftien  profiles  fall  below  "zero"  in  discussions  and  diagrams  which 
follow,  the  implication  is  that  the  particle  emission  or  radiation  upon 
which  the  measurement  is  based  is  not  detectable  above  background  effects 
(which  are  subtracted  before  analysis).  That  is,  the  concentration  of  the 
element  of  interest  is  below  the  detection  limit  for  the  method  used. 

Although  it  is  possible  to  measure  absolute  concentrations  by  making 
use  of  known  reaction  oross  sections,  detector  efficiencies,  and  geometric 
factors,  measurements  can  be  more  accurately  and  conveniently  made  by  first 
obtaining  a  depth  profile  of  relative  concentration  and  then  making  a  sec¬ 
ond  measurement  (with  only  the  target  changed)  an  a  specimen  with  a  known 
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number  of  target  atoms  (of  the  element  of  interest)  all  in  the  near-surf aoe 
region. 


The  raw  data  obtained  in  NBP  measurements  is  in  the  form  of  a  set  of 
partiole  or  photon  counts  expressed  as  a  funotion  of  incident  ion  energy. 
This  funotion  is  an  approximate  profile  (especially  for  narrow  resonances ) 
and  is  frequently  adequate  for  many  uses  (when  converted  to  concentration 
as  a  function  of  depth).  However,  one  of  the  purposes  of  the  work  upon 
which  this  report  is  based  was  to  establish  depth-concentration  profiles 
that  are  nearly  independent  of  background  effects ,  ion-beam-energy  distri¬ 
butions,  ion  energy-loss  straggling  and  resoncnce  widths. 

Following  this  introductory  section,  there  will  be  a  disoussion  of  the 
principles  of  NRP,  the  methods  of  compution,  a  discussion  of  existing  soft¬ 
ware,  the  conversion  of  relative  to  absolute  concentration  profiles,  and 
oonoluding  remarks  with  suggestions  for  future  development. 


; 
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II.  PRINCIPLES  OF  NRP 


This  discussion  of  the  principles  of  NRP  will  be  aided  by  an  illustra¬ 
tive  example,  a  measurement  (made  in  January  1979)  of  hydrogen  concentra¬ 
tion  in  an  SiOg-Si  specimen  (target)  by  means  of  the  reaction  ?Li(H,  V)8  Be. 
The  atoms  of  the  element  that  participates  in  the  reaction  upon  which  the 
measurement  is  based  will  be  called  the  target  atoms.  The  method  described 
need  not  be  confined  to  (p ,?)  reactions;  (p ,cQ  reactions  have  also  been 
used. 


A.  General  Features  of  the  Example.  Figure  1  shows  the  relative  gamma-ray 
yield  (, counts  per  minute  in  a  specified  energy  window)  from  the  7Li(H,'Y)®Be 
reaction  (circles  and  left  linear  vertical  scale)  and  the  absolute  hydrogen, 
concentration  (solid  histogram)  as  functions  of  the  7Li*  incident  energy 
relative  to  the  resananoe  energy  (3,050  keV)  for  the  target  sketched  in  the 
upper  part  of  the  diagram.  Note  that  the  concentration  scale  and  profile 
sure  truncated  at  the  lower  end. 

The  angstrom  depth  scales  are  fixed  by  taking  the  ion-energy  (above 
resonance)  to  be  numerically  equal  to  the  energy  loss  (in  keV)  at  the  depths 
indicated.  The  specific  loss  rates  (stopping  powers)  are  those  of  Si  and 
Si02  (presence  of  H  ignored)  and,  therefore,  the  angstrom  depth  scale 
factors  for  SiOg  and  Si  are  different.  The  specific  energy  loss  rates  are 
considered  to  be  constants  for  the  depths  spanned.  The  thickness  (depth) 
of  the  SiOg  layer  was  known  from  sample  fabrication  methods  but,  in  princ¬ 
iple,  could  have  been  measured  by  methods  described  in  this  report. 

The  difference  in  shape  between  the  profile  and  the  gamma-ray  yield 
curve  is  due  to  the  vacuum  beam-energy  distribution  (VBED),  the  cross-sec¬ 
tion  shape  and  the  beam-energy-los3  straggling.  By  careful  analysis,  the 
true  profile  can  be  determined  from  the  gamma-ray  yield  ourves  to  an  approx¬ 
imation  dependent  upon  the  statistical  merits  of  the  data. 

B.  Beam-Energy  Distributions.  The  VBED  was  determined  in  the  measurement  of 
the  example  by  a  magnetic  beam  analyzer  followed  by  an  electrostatic  beam 
analyzer  (EBA).  The  theoretical  VBED  is  triangular  in  shape  but  (partly 
because  of  slit-edge  soattering)can  be  well  represented  by  a  Guassian  with  its 
mean  energy  equal  to  the  Van  de  Graaff  terminal  voltage  (singly  charged  ions) 
and  its  full  width  at  half  maximum  (FHHM)  determined  by  the  EBA  aperture 
settings  (2).  In  the  example,  the  VBED  is  represented  by  a  Guassian  with  a 
FHHM  of  1.5  keV.  In  later  parts  of  the  report,  the  VBED  will  be  represented 
by  s(%n»®b)  where  is  the  mean  energy  and  g  dB^  is  the  probability 

that  a  beam  particle  will  be  in  energy  interval  dE^  at  E^.  The  method 
of  analysis  will  hold,  of  course,  for  a  VBED  of  any  shape.  For  convenience, 
the  VBED  may  be  normed  to  unit  area.  That  is,  the  probability  that  a 
projectile  will  be  in  some  energy  interval  of  the  representation  is  unity. 

As  the  beam  ions  enter  the  target,  they  lose  most  of  their  electrons; 
the  average  energy  is  degraded  (at  pertinent  average  energy  loss  rates); 
the  distributions  are  broadened  (increasingly  with  depth)  and  skewed  by 
statistical  fluctuations  in  energy  loss  (straggling). 


relative  GAMMA-RAV  yield 


-4  K  1400  A  -SLOj, 


7  +\ 

Lc  LONS 

1 

-* - SILICON - j 

(SAMPLE  HH0S5)  / 

I _ I  — 1 - 1 - -t- 

o  400  800  |200  1600 

DEPTH  IN  Si  (JO3*) 


\l  E.NE  RGrY  RELATIVE  TO  RESONANCE  £N  ERG  Y  CkeV) 
AMD  TARGET  DEPTH  in  ENERGY  LOSS  UNITS  (ke*0 

Fig.  1  -  Relative  garana-ray  yield  (circles)  and  H-atom  concentration 
(solid  line)  for  sketched  target 


We  may  represent  the  internal  beam-particle  energy  distribution  (EPED) 
at  depth  x  within  the  target  by  h(Eb,E,x),  hdE  being  the  probability  that 
a  particle  will  be  in  energy  interval  dE  at  E.  An  IPED  is  broader  than  the 
VBED  and  has  a  mean  energy  which  is  lower  by  the  average  energy  loss  of  the 
projectiles  over  the  paths  to  depth  x.  The  function  h(Bh»E,x)  can  be  cal¬ 
culated  with  the  help  of  a  straggling  function  f (E«B^«x)  where  fdE  is  the 
probability  that  a  particle  in  energy  interval  dEin  at  Ej_a  before  entering 
the  target  will  be  degraded  into  energy  interval  dE  at  E  as  it  is  slowed 
(mostly  by  electronic  stopping)  enroute  to  depth  x.  We  then  have 

h(Bb^,x)  5  f  g^.E^)  fCE.E^.x)  dEin  (1) 

'Bin 

where  the  integral  extends  over  the  VBED. 

C.  Gamma-ray  Yield  Functions.  If  the  incident  ion  energy  is  increased  in 
suitable  steps  from  below  resonance  energy  to  well  above  resonance  energy, 
and  a  photon  counter  (with  an  appropriate  energy  window)  intercepts  a  con¬ 
stant  fraction  of  the  emitted  gamma  rays,  a  relative  gamma-ray  yield  func¬ 
tion  such  as  that  plotted  (circles)  in  Fig.  1  will  be  obtained.  The  IfRP 
method  is  most  effective  when  there  is  a  strong,  narrow,  isolated  resonance 
for  capture  of  the  incident  atoms  into  a  compound  nucleus.  A  narrow  res¬ 
onance  provides  good  similarity  between  the  profile  and  the  yield  shape  at 
shallow  target  depths  where  straggling  is  small  and  provides  for  better 
ultimate  depth  resolution  than  does  a  broad  resonance.  It  is  assumed  for  this 
discussion  that  there  are  no  resonance®  associated  with  the  host  materials 
in  the  incident  energy  region  of  interest.  If  there  are,  the  associated 
yield  (as  well  as  any  nonresonant  yield  from  the  host  and  target  atoms) 
must  be  subtracted  before  futher  data  analysis. 

The  absolute  gamma-ray  yield  function  for  an  "impurity”  element  in  a 
homogeneous  host  of  thickness,  Tx»  can  be  represented  by 

rabs<Vrx>  S  h nef  /T*V(*)  MB^E.x)  ^(Er,E)  dE  dx  (2) 

Where  (  is  the  efficiency  and  solid-angle  faotor  of  the  detector;  If, 
is  the  number  of  incident  ions;  If ^<p(x)  is  the  concentration  (atoms/om^) 
function  for  the  element  being  profiled;  i^E^EfX)  is  the  IPED?  is 
the  total  number  of  atams/cm2  to  depth  T*  cr(E,Er)  is  the  cross-section 
function.  Sinoe  only  the  relative  yield  is  of  interest  in  the  description 
ahead,  the  subscript  abs  and  the  factors  %n  and  e  may  be  dropped.  The 
functions  under  the  integral  signs  may  be  normed  as  convenience  dictates. 
Absolute  profiles  are  derived  frcm  relative  profiles  by  means  of  data  frcm. 
a  calibration  target  of  known  composition;  this  will  be  discussed  in  Sec¬ 
tion  7.  Although  targets  may  be  composed  of  several  homogeneous  layers  of 
different  materials,  the  analysis  must  deal  with  each  homogeneous  layer 
separately  beoauee  of  its  particular  speoific  average  energy  loss  rate. 

Note  that  yield  is  a  funotion  of  target  thickness  Tx.  For  a  target  of 
moderate  thickness,  this  dependence  is  negligible  (for  a  narrow  resonance) 
except  near  the  end  of  the  target. 
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The  integrand  of  equation  (2)  can  be  viewed  aa  a  product  of  three  prob¬ 
ability  factors  (all  normed  to  unit  area  for  this  paragraph).  JL^x)  ia 
the  probability  that  there  ia  one  atom  per  unit  area  in  target  element  dx 
at  depth  x;  h(B^,B,x)  i3  the  probability  that  for  incident  energy  B«j  there 
will  be  a  projectile  at  depth  x  with  energy  E;  <7(3r,E)dE  i8  the  probability 
for  gamma-ray  emission  due  to  the  element  dE  at  E  of  the  cross-8eotion 
function  if  there  ia  one  target  atom  per  unit  area  in  element  dx  at  x. 

In  the  analysis  of  NRP  data,  it  is  oonvenient  to  express  target  depth 
in  terms  of  .d(x),  the  average  energy  loss  of  a  projectile  in  penetrating  to 
depth  x,  rather  than  in  terms  of  angstroms  or  am.  If,  in  diagrams  such  aa 
that  of  Pig.  1,  the  incident  energy  E^  a  E  is  associated  with  zero  depth 
in  the  target,  and  the  incident  energy  is  increased  from  B^  to  E^  +1  ,  the 
average  energy  of  the  IPED  at  depth  x  will  be  Er»  The  bombarding  energy 
(relative  to  resonance  energy)  is  numerically  equal  to  the  target  depth  in 
energy-loss  units. 

Equation  (3)  is  an  expression  for  relative  yield  with  the  depth  variable 
in  energy-loss  units. 

Y(Bb,T)  =  f  f  fT  H(I)  8(3^,^)  P(Ein,3,i)  o(Er,E)  dE  dEin  d 5  (3) 

•'Bjn-'E  b 

The  order  of  integration  is,  of  course,  immaterial  but  it  is  often  helpful 
to  integrate  with  respect  to  dB^j  first  to  give  the  analog  of  equation  (2). 

To  recapitulate  then,  we  see  that,  in  order  to  calculate  the  relative 
ganmia-ray  yield  function  Y(E^,T),  we  need  a  concentration  function  N(j), 
a  VBED,  a  straggling  function  P(Ein,E,3),  and  a  cross-section 

f inctiono(E,Br)._  In  the  example,  a  set  of  straggling  functions  (for  a 
discrete  set  of  3 )  is  calculated  from  the  theory  of  Vavilov  (3).  They 
could  have  been  calculated  from  other  theories  and  (for  economy  of  comp¬ 
utation  time)  should  be  calculated  at  different  depths  in  the  target  by 
different  methods;  for  example,  by  Vavilov’s  theory  for  shallow  depths 
(for  low-Z  targets  and  projeotiles  only)  and  by  Bohr  theory  (4)  for  great¬ 
er  depths.  Tifhile  the  correspondence  between  yield  and  profile  is  most 
obvious  for  narrow  ,  isolated  resonances,  there  is  no  reason  why  the 
cross-section  function  cannot  span  more  than  one  resonance.  The  detector 
energy  window  must  be  such  that  as  the  set  of  E^  is  traversed,  there  is 
no  yield  from  resonances  other  than  those  represented  by  the  cross-section 
function.  Recall  also  that  any  yield  from  the  host  must  be  considered  as 
background. 
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III.  PROFILE  DETERMINATION. 

The  relative  gamma-ray  yield  can  be  measured  and  can  be  calculated 

from  equation  (3)  if  the  profile  N(3)  is  known.  However,  the  objective 

is  to  measure  function  Y(Ey),T)  and  from  this  measurement  determine  N(3). 

This  may  be  done  by  guessing  at  a  profile  after  inspecting  the  measured 
yield  function,  calculating  Y(3^,T),  comparing  measured  and  calculated 
yields,  then  adjusting  N(5')  and  iterating  until  the  calculated  and  measur¬ 
ed  yield  functions  are  judged  to  be  in  reasonable  agreement.  Machine  cal¬ 
culations  and  plotting  facilitate  this  method.  For  profiles  of  the  shape 
of  that  of  Fig.l,  good  agreement  can  be  obtained  in  about  three  trials. 

To  carry  out  numerical  analyses,  the  entire  target  will  be  divided  into 

elemental  slices  perpendicular  to  the  depth  axis.  These  slices  will  be 
labelled  with  index  i  (i  *  1,2,3,»**»  )  with  the  surface  slice  labelled 
with  i  s  1.  The  slice  widths,  Wj_,  need  not  be  ’uniform.  This  targat  sub¬ 
division  will  be  independent  of  any  subdivision  into  homogeneous  layers  of 
different  stopping  power  (layered  targets).  Let  us  first  consider  the 
target  to  be  homogeneous  and  later  consider  steps  needed  to  treat  layered 
targets . 

In  applications  of  the  NRP  method,  the  need  frequently  arises  to 
determine  the  profiles  of  an  "impurity”  element  in  many  identical  hosts. 

For  example,  in  solid-state  device  development,  it  may  be  useful  to  know 
the  hydrogen  profiles  in  the  sample  configuration  of  Fig.  1  for  many  ident¬ 
ically  prepared  samples  which  have  been  subjected  to  a  varied  sot  of  process¬ 
ing  environments  (annealing  etc.).  For  this  reason  and  as  a  basis  for  the 
possible  development  of  a  more  sophisticated  analysis  method,  it  is  use¬ 
ful  no  construct  from  equation  (3)  a  matrix,  M,  characteristic  of  the  reac¬ 
tion  parameters,  the  mode  of  target  subdivision,  but  independent  of  the 
profile.  This  matrix  together  with  a  vector,  P,  characteristic  of  the  profile 
and  the  mode  of  target  subdivision  will  permit  calculation  of  the  yield  func¬ 
tion  according  to 


Y(Eb,T)  =  MP  (4) 

where  Y(E{j,T)  is  a  set  of  ordinates  of  the  calculated  yield  function  at 
a  set  of  absoissas,  -  By,  the  latter  set  being  numerically  equal  to  an 
energy-loss  set  corresponding  to  various  depths  into  the  target. 


'Biis  suggests  a  method  of  profile  determination  without  the  need  for 
iteration.  If  M  can  be  inverted  to  M"l,  then 


P  =  IT1!  . 


(5) 


Some  exploratory  work  has  been  done  toward  this  goal  but  an  algorithm 
for  the  inversion  which  will  produce  nonsingular  inverse  matrices  with 
acceptable  errors  has  not  yet  been  devised. 
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The  ith  column  of  the  matrix,  M,  is  obtained  by  numerically  integrat¬ 
ing  equation  (3)  with  respect  to  dE.  and  dE.  The  row  indices,  m,  are 
those  of  the  set  E^.  The  column  indices,  i,  are  those  of  the  slices. 

The  ith  element  of  P  is  N(-A)w^.  The  yield  associated  with  (B^)m  is  then 


LsTb  .P.  .  (6) 

a  £  m,i  i  v  ' 

For  targets  composed  of  more  than  one  homogeneous  layer,  (five  is 
the  maximum  in  the  present  software)  the  yield  for  a  particular  value  of 
(Eb)jn  obtained  by  summing  the  calculated  yields  for  each  layer,  com¬ 
puted  separately  with  the  IPSD  for  the  first  interface  substituted  for 
g(B^,Ein )  to  get  the  contribution  for  the  second  dayer  and  so  on.  That 
is,  the  projectile  energy  distribution  to  be  used  for  btaining  the  yield 
from  the  nth  layer  (n  s  2,3,«* **)  is  the  IPED  at  its  most  shallow  boundary. 

The  target  subdivision  scheme  should  be  such  that  the  number  of  sub¬ 
divisions  is  kept  to  a  minimum  with  the  proviso  that  a  finer  subdivision 
will  not  alter  the  calculated  yield  function  significantly.  In  some  cases, 
it  may  be  useful  to  choose  the  width  of  the  first  slice  to  be  as  small  as 
ten  angstroms  if  the  incident  energy  distribution  and  the  resonance  are 
very  narrow  and  precise  information  very  near  the  surface  of  the  target  is 
wanted.  The  slice  widths,  w^,  nay  increase  approximately  as  the  square 
root  of  the  depth  (without  losing  obtainable  details  of  the  calculated 
yield  function)  since  straggling  effects  increase  in  this  manner  and  they 
tend  to  smear  sharp  details  in  the  parts  of  the  measured  yield  curve  assoc¬ 
iated  with  the  interior  slices  of  the  target. 

Once  the  target  subdivision  scheme  has  been  chosen,  a  set  of  func¬ 
tions  Fi(<l,3i)  must  be  calculated,  one  for  each  slice;  \  is  the  exact  energy 
loss  of  a  projectile  enroute  from  the  surface  of  a  target  layer  to  the 
center  of  the  ith  slice;  3 i  is  the  average  energy  loss  in  keV  from  the 
layer  surface  to  the  center  of  the  ith  3lice.  The  present  software  calcul¬ 
ates  the  set  by  means  of  Vavilov's  theory  (3)  which  has  limited  appli¬ 
cability  as  will  be  discussed  later. 

In  our  example,  g(Ej5,Ej_n)  is  represented  by  a  Gaussian  function  and 
o(E,Er)  by  a  Breit-Wigner  sincle-resonance  function. 

Many  NRP  calculations  and  plots  have  been  made  with  the  NPL  Texas 
Instruments  Advanced  Scientific  Computer  (ASC7)  via  a  keyboard  terminal. 

Data  and  program  files  associated  with  the  illustrative  example  are  in 
the  ASC7  catalog  at  node  USERCAT/D66/B^/DUlIMl/NRP,,,  which  will  be  replac¬ 
ed  with  the  JSL**  synonym,  P5,  in  the  following  discussion.  The  source 
language  is  Fortran. 

Some  detail  concerning  the  software  will  appear  in  the  next  section 
but  what  follows  immediately  will  simply  show  how  to  use  files  and  programs 
pertinent  to  the  example. 


*  The  symbol  jS  represents  sero. 

**  ASC7  job  specification  language. 
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The  files  needed  to  run  the  calculation  of  the  example  are  sons  of 
the  catalog  node  P5;  that  is,  e.g.,  P5/SL9YMAT  is  the  catalog  path  to  the 
source-language  file  of  the  program  Yl'AT  that  calculates  the  matrix. 

The  object-language  file  and  the  load  module  for  YMAT  are  the  sons  0L9YMAT 
and  LM9YMA.T  respectively.  The  input  data  file  (son)  is  IN2TGT7  and  the 
output  file  containing  the  matrix  is  MA2TGT7.  Usually  the  scheme  has  been 
to  label  NHP  files  with  the  first  two  characters  indicating  the  kind  of 
file,  with  the  next  character  as  a  "serial  number"  indicating  a  particular 
version,  and  with  the  last  four  characters  indicating  that  the  file  is  assoc¬ 
iated  with  a  particular  program  or  target. 

The  computer  run  of  the  program  YMAT  may  be  made  with  the  help  of 
the  file  at  P5/HLP7,  suitably  modified.  This  is  a  JSL  file  and  may  be 
run  from  a  keyboard  terminal,  piecewise  if  desired.  There  is  provision 
for  a  file  named  FT59F0j?(l  which  contains  diagnostic  plots  ihich  may  have 
been  initiated  in  the  input  file,  IN2TGT7,  at  the  user's  option. 

Before  making  use  of  the  matrix  of  our  example,  let  us  turn  to  the 
processing  of  the  experimental  data  and  the  cataloging  of  the  results 
in  a  form  that  will  be  suitable  for  comparison  with  calculated  yield. 

The  results  of  measurements  have  been  formatted  and  filed  at  P5/IXHE0S5. 

A  JSL  program  similar  to  that  at  P5/HLP5  facilitates  the  processing  of 
experimental  data  via  program  DAIA5  with  source-language,  object-language, 
and  load-module  files  at  P5/SL1DATA5,  P5/0L1DATA5,  and  P5/LM1DATA5  res¬ 
pectively.  Program  DATA5  provides  for  a  plot  of  the  experimental  data 
that  will  be  needed  for  the  later  conversion  of  a  relative  profile  to  an 
absolute  profile.  The  output  file  of  DATA5  is  cataloged  at  P5/0X1HK0S5. 

The  next  step  is  to  prepare  a  file  containing  a  trial  profile,  merge 
it  with  the  experimental  data  file,  0X1HH0S5,  and  with  the  file  MA2TGT7. 

This  merged  file  for  the  example  is  at  P5/IMHH0S5. 

This  merged  file  is  the  input  file  for  program  YMULT  that  performs 
the  matrix  multiplication  and  produces  a  plot  of  the  essential  information 
of  Fig.  1.  The  source,  object,  and  load-module  files  for  YMULT  are  at 
P5/SL8YMULT,  P5/0L8YMJLT ,  and  P5/LM8YMULT  respectively.  The  JSL  list  at 
P5/HLP3  or  one  similar  to  it  oan  be  of  help  in  running  YVULT. 

If  the  plotted  calculated  yield  does  not  compare  favorably  with  the 
plotted  measured  yield,  the  profile  of  file  IMHE0S5  can  be  altered  and 
another  YMULT  run  made. 
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17.  EXISTING  S0FT7TAHE 


Fortran  programs  for  carrying  out  the  kind  of  analysis  described  in 
preceding  sections  of  this  report  reside  in  the  catalog  of  the  HKL  Texas 
Instruments  ASC7  computer.  Program  listings  will  not  be  given.  Instead, 
catalog  path  names  to  pertinent  files  will  be  given  and  programs  will  be 
discussed  with  the  help  of  the  example  of  Fig.  1. 

The  source-language  programs  are  annotated  so  that,  hopefully,  they 
can  be  understood  and  used  with  a  reasonable  effort.  These  programs  were 
written  with  3mall  effort  toward  sophistication  but  rather  with  the  goal 
of  creating  and  using  the  software  in  applications  of  some  urgency.  There¬ 
fore,  they  are  not  in  the  exceedingly  rare  category  of  programs  that  cannot 
be  modified  to  better  serve  the  user. 

The  catalog  path  names  for  files  pertinent  to  this  report  will  include 
the  JSL  synonym,  P5  *  USERCAT/D66/B^?/DIJITNK1/NRP.  A  catalog  3tatus  print¬ 
out  for  the  node  associated  with  this  path  appears  in  the  Appendix.  The 
JSL  listings  at  P5/HLP1,  P5/HLP3,  etc.  may  be  of  some  help  in  running  and 
in  modifying  the  programs. 

A.  Initial  Data  Treatment.  The  data  obtained  from  measurements  is  enter¬ 
ed  into  a  file  to  be  read  by  program  DATA5.  The  names  of  the  quantities, 
their  meanings,  and  formats  oan  be  obtained  from  the  source-language  list¬ 
ing  at  P5/SL1DATA5  and  the  input  file  for  the  example  at  P5/IXHH0S5. 

Program  DATA5  first  reads  and  prints  the  input  file  which  (for  purpo¬ 
ses  of  this  discussion)  will  be  said  to  consist  of  card  images  named: 
l,2,3,4a,4b,*»*« ,  and  5.  On  the  first  ran,  set  NBHREJ  s  p  in  image  1  and 
omit  image  2.  If  the  results  of  the  first  run  show  that  there  are  N  ob¬ 
viously  wild  data  points  among  those  represented  in  images  4a,  4b,  •  •», 
set  NBHHEJ  “  Nj  enter  the  ordinal  numbers  of  the  points  to  be  rejected  in 
image  2  and  make  a  second  run  with  the  amended  file.  If  NBR5EJ  p  p,  sub¬ 
routine  CHANGE  will  be  called  to  make  the  deletion.  Experience  indicates 
that  this  feature  of  the  program  is  almost  never  used. 

Card  image  3  (second  image  in  file  IXHH0S5)  may  contain  up  to  SO 
Hollerith  comment  characters.  In  the  input  file  for  the  example,  "NO 
TILT"  means  that  the  beam  was  striking  the  target  90  degrees  to  the  plane 
of  the  target;  (some  smaller  angle  could  be  used  to  enhance  depth  resolu¬ 
tion).  The  notation  "NV7CT  =  79.0"  indicates  that  the  controlled  temperature 
of  the  SBA  (northwest  corner  temperature)  was  79.0  degrees  Fahrenheit. 

Card  image  4a  contains  data  for  the  first  three  measurement  points, 
the  data  for  each  point  consisting  of  an  ordinal  number  (NKUN),  the  time 
of  the  run  in  minutes,  the  number  of  yield  counts,  and  the  potentiometer 
setting  (of  the  SBA).  These  data  sire  entered  in  order  of  increasing  pot¬ 
entiometer  (beam  energy)  settings.  Card  image  5  contains  six  numbers. 

The  first  is  the  potentiometer  constant  (POTCST  *  POTSET  *  1000  gives  the 
beam  energy  in  keV). 
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BGlCST  is  the  number  of  counts  per  minute  obtained  in  a  background  run 
v/ith  the  beam  hitting  a  stopper  (distant  from  the  target).  SDBG1  is  cal¬ 
culated  by  the  user  according  to  the  formula  SQRT(BGlCST/(number  of  min¬ 
utes  for  background  run)).  This  quantity  is  the  standard  deviation  of 
BGlCST. 

The  third  number  in  card  image  5  permits  the  average  yield  for  the 
first  NBG2  points  to  be  subtracted  from  all  points;  NBG2  s  0  means  no 
subtraction.  This  constant  facilitates  the  subtraction  of  a  pre-reson¬ 
ance  beam-dependent  background  under  the  following  circumstances.  If  a 
single*  narrow,  isolated  resonance  is  the  basis  for  the  measurement  and 
the  bombarding  energy  for  the  first  point  is  far  below  the  resonance  energy, 
the  yield  for  the  first  few  points  will  be  essentially  that  due  to  the  host 
and  the  tails  of  any  distant  resonances;  this  background  should  be  subtract¬ 
ed  from  all  points.  For  example,  if  the  yield  from  the  first  three  data 
points  is  small  and  nearly  equal,  set  1JBG2  Z  3.  No  subtraction  of  this 
kind  was  made  for  the  analysis  of  the  example. 

The  last  three  numbers  in  card  image  5  have  meanings  us  follows. 

IJPLOT  =  1 , ($)  causes  (prevents)  the  calling  of  subroutine  GRAPH3  which 
creates  a  plot  file  (FT59F(2^1)  which  can  be  used  to  actuate  a  Calcomp 
plotter  to  produce  a  plot  of  the  experimental  data.  This  plot  will  show 
the  raw  yield  counts  (squares),  the  yield  after  time-dependent  background 
subtraction  (crosses),  the  yield  after  target-dependent  background  sub¬ 
traction  (triangles),  and  vertical  lines  indicating  the  standard  deviation. 
NPNCH  s  1,(^0  causes  (prevents)  the  punching  of  the  output  file  onto  cards 
(file  0X1HH0S5  for  the  example).  fJFILB  ■  1,(0)  causes  (prevents)  the  print¬ 
ing  of  the  output  file  when  DATA5  is  being  run. 

The  plot  of  experimental  data  i3  useful  (as  will  be  explained  in  Sec¬ 
tion  V)  for  converting  relative  concentration  profiles  to  absolute  concen¬ 
tration  profiles. 

A  more  accurate  method  of  determining  the  beam-dependent  background 
consists  of  obtaining  a  measured  yield  curve  for  the  host  devoid  of  the 
impurity  (element  of  interest)  for  an  identical  set  of  bombarding  energies 
and  making  a  point-by-point  subtraction  to  provide  input  data  for  program 
DATA5.  This  is  especially  feasible  when  the  impurity  is  an  implanted  one. 

B.  Matrix  Computation.  The  matrix,  M,  of  equation  (4)  is  computed  by 
program  YMAT  with  source-language  file  at  P5/SL9YMAT,  object  file  at 
P5/0L9YMAT,  and  load  module  at  P5/LM9YMAI.  This  matrix  is  dependent  upon 
beam  parameters,  resonance  parameters,  and  all  target  parameters  except 
the  concentration  function  (profile)  of  the  element  of  interest.  It  is 
obtained  by  integrating  equation  (3)  with  respect  to  E  and  Ein.  The  inte¬ 
gration  with  respect  to  A  is  replaced  by  matrix  multiplicationcarried  out 
by  another  program  named  YMXT. 

The  YMAT  source-language  listing  contains  many  comments  and  a  card- 
image  table  for  the  input  file.  A  convenient  way  to  create  a  new  input  file 
for  YMAT  is  to  display  an  existing  input  file  (for  example,  that  for  our 
illustrative  example  at  P5/IN2TGT7)  on  a  keyboard  terminal  and  alter  it  as 
needed. 
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The  source-language  listing  may  be  used  to  obtain  an  understanding  of 
program  YLIAT.  Additional  comments  which  will  be  of  help  follow. 

The  statement  labels  in  the  Fortran  listing  are  not  always  in  ascend¬ 
ing  order;  this  is  due  in  part  to  the  use  of  blocks  of  statements  from 
earlier  programs. 

Program  YMAT  first  reads  and  prints  the  input  file  (IN2TGT7  for  our 
example). 

The  next  action  is  to  compute  the  parameters  needed  for  the  division 
of  the  target  into  slices;  later,  an  energy-loss  distribution  function 
(ELD)  will  be  associated  with  each  slice.  The  average  loss  for  a  particu¬ 
lar  ELD  will  be  that  to  the  center  of  the  associated  slice.  The  slices  at 
the  surface  of  the  target  will  have  thioknesses  a  few  times  the  FWHM  of 
the  VBED.  The  interior  slices  will  increase  in  thickness  with  depth  into 
the  target  at  a  rate  roughly  proportional  to  the  increase  in  energy-los3 
straggling  with  depth  in  order  to  minimize  the  number  of  slices  without 
sacrificing  depth  resolution.  That  is,  the  interior  resolution  is  limited 
by  the  of  the  IPED  which  increases  with  depth  due  to  straggling. 

The  widths  of  the  ELD  increase  roughly  as  the  square  root  of  depth  of 
penetration.  This  subdivision  into  slices  spans  the  various  homogeneous 
layers  that  nay  exist. 

After  reading  and  printing  the  input  file,  program  YJJAT  calls  sub¬ 
routine  LAMGR  (layer  manager)  to  select  the  current  (initially  the  first) 
homogeneous  layer. 

The  next  task  is  to  compute  parameters  needed  in  the  calculation  of 
the  KT.D  set.  The  calculation  of  these  functions  is  based  upon  a  theory 
due  to  P.  V.  Vavilov  (3)  and  a  modification  of  a  program  kindly  supplied  by 
S.  M.  Seltzer  and  M.  Berger  (5).  A  discussion  of  the  suitability  of  this 
theory  for  the  purpose  at  hand  will  be  found  in  the  final  section  of  this 
report.  The  loss  axis  for  the  EUD  set  must  span  a  domain  great  enough  to 
accommodate  the  greatest  loss  suffered  in  the  target  by  any  beam  particle; 
in  application,  this  means  that  the  loss  domain  must  extend  from  zero  to 
some  value  somewhat  greater  than  the  target  thickness.  For  purposes  of 
ELD  calculation,  the  loss  axis  is  divided  into  uniform  intervals  of  width 
DE  with  ordinal  numbers,  K,  (sometimes  J5);  the  interval  of  greatest  loss 
is  designated  by  K  s  1;  the  ordinal  numbers  increase  with  decreasing  loss. 

To  save  computer  costs,  the  ordinates  of  the  FT.n  are  not  calculated 
in  detail  for  every  loss  interval,  DE,  except  in  the  low-loss  region. 

Rather,  the  loss  axis  is  divided  into  blocks  and  the  ordinates  for  sane 
intervals  are  supplied  by  interpolation  as  explained  in  the  Fortran  list¬ 
ing.  A  maximum  of  five  blocks  is  provided  for;  the  detailed  calculations 
are  leS3  danse  as  the  block  numbers  increase. 

The  independent  variable  in  the  Vavilov  functions  is  not  energy  loss 
but  is  a  quantity  named  lambda  which  is  a  function  of  the  exaot  energy  loss 
and  the  average  energy  loss  of  interest.  These  funotions  must  be  converted 
into  functions  of  aotual  energy  loss  before  use  in  this  analysis. 


In  the  midst  of  the  ELD  parameter  calculations,  the  target  slice 
widths  are  calculated  by  calling  subroutine  TARG1  which  may,  in  turn, 
call  TAXG2  and  TARG3.  Upon  completion  of  these  tasks,  many  of  the  quan¬ 
tities  calculated  up  to  this  point  by  the  main  program  (YMAT)  are  printed. 

A  set  of  bombarding  energies,  E2(M),  is  needed  for  the  calculation 
of  the  matrix  elements;  the  rows  of  the  matrix,  M,  are  labelled  by  the 
integer  variable,  M;  this  set  need  not  intersect  the  set  of  bombarding 
energies  used  in  the  measurement.  The  YMAT  input  file  contains  parameters 
which  permit  the  set  EB(M)  to  have  intervals  of  three  different  sizes  in 
order  to  accommodate  the  sharp  rise  in  yield  at  the  target  surface  assoc¬ 
iated  with  narrow  resonances.  The  set  of  bombarding  energies,  the  resonance 
function,  and  the  VBED  are  calc  dated  by  calling  subroutine  SETS6. 

Subroutine  SETS6  calculates  a  Gaussian  VBED  using  the  FVHM  provided 
in  the  input  file.  This  function  is  arbitrarily  cut  off  at  1.5  FWHM 
since  the  exact  shape  of  the  VBED  is  not  known  (theoretically,  it  is  a 
triangle  modified  by  slit-edge  scattering)(2).  The  cross-section  function 
is  calculated  by  means  of  the  Breit-YTigner  single-resonance  formula  with 
the  skirts  cut  off  when  the  magnitude  of  the  abscissa  is  NGAK  FVHII,  where 
NGAM  is  from  the  input  file.  The  density  of  ordinates  calculated  is  decreas¬ 
ed  in  the  regions  remote  from  the  peak  of  the  resonance  function  to  save 
computation  time.  To  encompass  all  but  a  few  tenths  percent  of  the  area 
under  the  curve,  the  abscissas  must  be  extended  to  about  100  F7I3M  on  either 
side  of  the  peak.  In  the  present  calculations,  cutoff  can  be  made  much  soon¬ 
er  without  significant  effeot. 

The  beam  energy  distribution  and  the  resonance  function  can  be  plot¬ 
ted  for  diagnostic  purposes  at  the  option  of  the  user.  Subroutine  GRAPHC 
is  called  for  this  purpose  if  parameters  are  set  in  the  input  file  as 
explained  in  the  listing. 

Program  YMAT  is  now  ready  to  calculate  (by  columns)  the  elements  of 
the  matrix  M.  The  column  index  is  I  (vdiich  also  designates  a  target 
slice  and  the  ELD  associated  with  that  slice).  The  calculation  begins  with 
the  lowest  I-value  for  the  current  homogeneous  target  layer  and  with  the 
ELD  of  greatest  width  (smallest  K). 

Subroutine  VAV7  is  called  to  compute  the  ELD  for  the  Ith  slice  and, 
in  turn,  calls  subroutine  T3IGEX  to  supply  oertain  needed  functions. 
Parameters  in  the  input  file  (see  curd-image  table  in  the  listing)  permit 
considerable  diagnostic  printout.  Provision  may  be  made  (in  the  input 
file)  to  plot  (via  subroutine  GHAPH6)  selected  ELD  in  order  to  see  if  they 
have  reasonable  shapes,  smoothness,  and  cutoff. 

Next,  subroutine  IPED2  is  called  to  provide  an  IPED  shape  for  the 
•  current  target  layer.  This  is  the  VBED  shape  for  the  surface  layer;  for 

interior  layers,  this  shape  must  be  modified  to  take  into  acoount  the 
average  energy  loss  and  straggling  in  earlier  layers.  This  shape  is  then 
!  used  in  conjunction  with  the  set  of  bombarding  energies,  EB(li)  to  provide 

>  IPED  for  the  integration  with  respect  to  E  in  equation(2).  Note  that  the 

IPED  corresponds  to  the  h(Bb,E,x)  of  equation  (1). 
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Subroutine  TINT  is  called,  to  perform  the  integration  over  E  of  equa¬ 
tion  (3)  to  produce  the  matrix  element  for  current  M  and  I.  After  this 
is  done  for  all  M  (all  elements  in  the  Ith  column),  1  is  advanced  and  the 
3ubnatrix  for  the  current  layer  is  completed. 

The  submatrices  for  other  target  layers  are  computed  and  the  yield 
for  a  uniform  (unit  height)  concentration  function  is  computed  for  diag¬ 
nostic  purposes  and  may  be  plotted  at  the  option  of  the  user. 

The  YMAT  output  file,  -which  contains  the  matrix  and  other  information 
for  transmission  to  program  YMULT,  is  written  to  a  file  and  cataloged. 

That  of  our  example  is  at  P5/AlA2TGT7. 

This  program  may  be  run  on  the  ASC7  with  a  JSL  listing  similar  to 
that  at  P5/hLP7. 

C.  Yield  Computation.  The  program  YMULT  computes  the  yield  function  by 
multiplying  a  trial  profile  function  by  the  matrix,  M,  which  was  generated 
by  program  YMAT. 

The  YMULT  input  file  for  the  illustrative  example  is  cataloged  at 
P5/IMHH0S5.  The  first  two  card  images  contain  constants  that  control  the 
reading,  printing,  and  plotting  done  by  YMULT.  These  are  followed  by 
three  merged  files;  the  first  is  the  output  file  from  DATA5  which,  for  our 
example,  is  cataloged  at  P5/0X1HH0S5.  The  next  is  a  trial  profile  which, 
for  convenience,  has  a  maximum  value  of  unity.  The  normalization  can  be 
chosen  at  will  because  we  are  calculating  a  relative  profile.  Next  is  the 
output  file  of  the  program  YMAT  which,  for  our  example,  is  cataloged  at 
P5/ilA2TGT7.  This  YMULT  input  file  can  be  conveniently  formed  by  modifying 
an  existing  file  and  merging  files  at  a  keyboard  terminal. 

Program  YMULT  can  be  run  on  the  ASC7  with  a  JSL  listing  similar  to 
that  cataloged  at  P5/HLP5 .  The  source-language  listing,  object  listing, 
and  load  module  are  at  P5/SL8YMULT,  P5/0L8YMULT  respectively. 

The  following  remarks  together  with  the  Fortran  listing  will  be  of 
help  in  understanding  and  in  the  use  of  program  YMULT. 

There  is  provision  in  YMULT  to  read  and  prooess  more  than  one  trial 
profile  funotion  but  this  feature  has  been  seldom  used. 

The  actual  matrix  multiplication  is  accomplished  by  calling  subrou¬ 
tine  MULT.  The  rows  of  the  matrix  computed  by  YMAT  are  multiplied  into 
a  column  vector  with  each  element  composed  of  the  product  of  a  concen¬ 
tration  value  associated  with  the  Ith  target  slice  and  the  width  (in 
energy-loss  units)  of  the  Ith  slice  in  order  to  produce  a  yield  value 
for  each  bombarding  energy,  EB(M).  This  is  done  for  each  target  layer  and 
the  layer  yields  sunned  in  an  appropriate  manner. 

At  the  option  of  the  user,  YMULT  will  call  subroutine  G3AFHD  to 
create  a  file  for  the  Calcomp  plotter,  thus  providing  for  plots  of  the 
kind  shown  in  Fig.  1  -  except  that,  in  addition  to  the  profile  (whioh 
will  be  plotted  with  +  signs)  and  measured  gamma-ray  yield  (whioh  will 
be  plotted  with  octagons),  the  calculated  yield  will  be  plotted  (with 
triangle  symbols). 


If  inspection  of  this  plot  shows  that  the  agreement  between  the  meas¬ 
ured  and  calculated  yield  is  satisfactory  in  the  user's  judgement,  the 
trial  profile  plotted  can  be  considered  a  useful  approximation  to  the 
actual  profile.  If  not,  the  trial  profile  nay  be  altered  and  another  run 
of  YFJLT  made. 

To  facilitate  comparisons  of  the  measured  and  calculated  yields, 
the  experimental  yield  is  scaled  for  plotting  so  that  the  areas  under 
both  ourves  out  to  an  abscissa  EBNBM  (a  parameter  inserted  by  the  user 
into  the  IMCrLT  input  file)  are  equal. 

All  plotting  programs  described  in  this  report  were  written  for 
the  CDC  3800  computer  and,  therefore,  in  the  linking  run  (3ee  JSL  listing 
cataloged  at  P5/HLP4 )  the  module  at  USESCAT/1.IB/SUP/0BJECT  must  be  assign¬ 
ed  to  the  job  to  permit  the  ASC7  to  use  t  ie  CDC  3800  code. 


V.  C0NVS2S Id'  OF  RELATIVE  PROFILES  ?C  ABSOLUTE  PROFILES. 


Relative  1,'RF  concentration  profiles  may  be  converted  to  absolute  con¬ 
centration  profiles  by  means  of  the  absolute  interaction  cross-section 
function,  geometric  factors,  detector  efficiency,  and  other  parameters 
associated  with  a  measurement.  This  may  also  be  done  with  the  help  of 
measured  yield  functions  from  calibration  targets  of  known  target-atom 
concentrations.  The  latter  method  is  probably  more  convenient  and  accurate 
than  the  former  and  will  be  discussed  for  thin  and  thick  targets. 

A.  Thin  Targets.  Targets  that  produce  yield  curves  that  (for  narrow 
resonances)  rise  from  zero  yield  (with  background  subtracted)  at  low  bombard¬ 
ing  energies  and  fall  to  zero  again  at  higher  (attainable)  bombarding 
energies  will  be  termed  thin  targets.  A  yield  function  associated  with 
a  thin  target  will  be  referred  to  as  an  entire  function.  Consideration 
of  equation  (2)  shows  that  the  area  under  a  yield  curve  that  is  an  entire 
function  is  proportional  to  the  total  number  of  target  atoms/cm2,  Up. 

The  constant  of  proportionality  may  be  found  from  the  area  under  a  measured 
yield  curve  for  a  thin  calibration  target  with  a  different  but  known  Up. 

An  example  of  such  a  calibration  targot  is  a  metal  or  semiconductor  host 
implanted  with  a  known  number  of  target  atoms/cm2,  all  in  the  first  few 
thousand  angstroms  below  a  well  defined  surface.  Convenient  units  for 
such  constants  are 


(No.  of  target  atoms )/cm2  . 

(keV  counts )/microooulomb 

The  area  under  the  relative  profile  must  be  equal  to  N™  for  an  un¬ 
known  thin  target.  This  allows  a  scale  (atoms/cm3)  to  be  placed  upon  the 
relative  profile  showing  the  absolute  concentration  for  the  unknown  target 
as  a  function  of  depth, 

B.  Thick  Targets.  The  yield  for  a  thick  tartet  in  not  an  entire  function 
of  the  bombarding  energy  and  the  analysis  for  thin  targets  does  not  pertain. 
That  is,  the  yield  does  not  fall  to  zero  near  the  end  of  the  attainable 
bombarding  energy  span.  Procedures  for  establishing  concentration  scales 
by  means  of  thick  calibration  targets  based  upon  careful  considerations 
(not  recorded  here)  follow. 

The  measured  yield  curve  for  a  thick  calibration  target  of  known 
uniform  target-atcm  concentration  will,  after  an  initial  rise  from  zero 
(narrow  resonance  assumed),  achieve  a  constant  value  and  remain  there  for  the 
remainder  of  the  bombarding  energy  span.  Some  selected  area  entirely  under 
the  constant  part  of  the  ourve  (keV  counts)  can  be  used  to  determine  a 
constant  which  may  be  used  with  thin  targets  of  unknown  concentration  as 
previously  described,  provided  that  the  (No.  of  target  atoms )/cm2  is  that 
associated  with  the  segnent  of  the  target  spanned  by  the  selected  area. 


If  the  target  of  unknown  concentration  is  thick,  a  convenient  value 
of  T  must  be  selected;  the  associated  calculated  yield  curve  must  have 
this  value  of  T;  the  measured  yield  curve  must  be  converted  to  an  entire 
function  by  means  of  T  and  ths  calculated  yield  curve  associated  with  the 
relative  profile.  The  concentration  scale  may  then  be  determined  as  before. 
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VI.  CONCLUDING  H3I1ASSS 


This  section  contains  a  few  general  comments  concerning  nuclear 
resonance  profiling  computational  methods  which  may  be  useful  in  possible 
future  developments. 

A.  Software  Architecture.  As  stated  previously,  the  existing  software 
was  written  with  the'  object  of  producing  working  programs  for  urgent 
applications  and  without  great  concern  for  the  user's  convenience  or 
computational  costs.  As  the  programs  were  used,  it  became  evident  that  there 
is  much  room  for  improvement  -  especially  with  respect  to  the  transparency 
and  ease  of  application  for  the  user.  Computational  costs  can  be  reduced 
by  restructuring  software.  Attention  should  be  given  to  creating  modules  that 
can  be  run  conveniently  on  computers  smaller  than  the  ASC7. 

One  fertile  area  for  improvement  concerns  the  computation  of  energy- 
loss  distribution  functions  (ELD).  Presently,  a  slight  change  in  target 
configuration  requires  essentially  the  sane  set  to  be  recomputed.  In 
addition,  all  vi.U  are  now  computed  from  Vavilov's  theory  and  interpolation 
methods  were  introduced  to  keep  computation  times  within  reasonable  limits. 

A  great  improvement  could  be  made  by  computing  sets  of  ELD  and  cat¬ 
aloging  them  in  computer  files  in  a  form  suitable  for  use  by  a  program 
like  YMAT.  Only  a  few  files  of  this  kind  would  be  necessary  for  measure¬ 
ments  made  in  connection  with  solid-state  device  development,  where  only 
a  few  host  elements  and  compounds  are  of  interest  such  as  Si  and  Si02  for 
example.  These  files  might  also  be  useful  for  other  work  involving  energy- 
loss  straggling. 

Various  computational  methods  could  be  used  in  creating  a  particular 
set.  Considerable  care  and  time  could  be  spent  on  a  set  that  is  to  be 
used  many  times.  These  sets  could  be  used  for  other  purposed  such  as, 
for  example,  in  particle  scattering  analysis.  As  new  theories  and  improved 
data  appear,  improved  sets  could  replace  existing  sets  without  need  to 
modify  using  programs. 

Vavilov's  theory  agrees  very  well  with  experiment  for  penetrations 
of  a  few  thousand  angstroms  for  elements  with  atomic  numbers  of  15  and 
below;  agreement  with  experiment  is  not  as  good  for  elements  of  much 
higher  atomic  number.  Some  other  theory  or  an  approximation  of  this 
theory  should  be  used  for  penetrations  of  more  than  a  few  thousand  ang¬ 
stroms,  e.g.,  Bohr's  theory.  Vavilov's  theory  relies  on  Bethe's  formula 
for  average  energy  loss.  The  algorithm  in  IMAT  avoids  this  dependence, 
and,  instead,  is  based  upon  measured  average  energy  losses. 

There  are  many  theories  relating  to  energy  loss  and  energy^loss 
straggling;  however,  there  are  few  measurements  available  to  test  strag¬ 
gling  theories.  These  measurements  are  difficult  and  are  badly  needed  for 
work  of  the  kind  described  is  this  report. 


B.  Computers  and  Peripheral  Apparatus.  Keyboard  terminals  located  in  the 
space  of  the  user  provide  a  significant  advantage  over  the  batch-processing 
mode  at  a  centrally  located  computer  for  calculations  of  the  kind  discuss¬ 
ed  in  this  report.  Keyboard-printer  terminals  can  be  augmented  with  video 
terminals,  plotters,  and  line  printers  dedicated  to  a  particular  set  of 
problems  thus  further  reducing  travel  and  waiting  times  associated  with 
batch-processing  modes  at  a  central  computational  facility. 

A  further  step  in  the  decentralisation  of  facilities  which  can  increase 
the  speed  and/or  ease  of  analyses  is  to  provide  the  user  with  a  small  com¬ 
puter  with  appropriate  peripherals  dedicated  to  a  particular  set  of  problems. 
That  support  for  such  arrangements  is  prevalent  can  be  inferred  from  the 
many  specialized  computers  dispersed  throughout  some  scientific  laborator¬ 
ies.  Of  course,  many  of  these  units  are  used  for  on-line  operations  but 
they  can  also  be  used  to  great  advantage  for  computational  purposes,  espec¬ 
ially  in  conjunction  with  a  large,  fast  machine  with  large  disk  storage 
capacity. 

This  mode  of  computation  has  been  abetted  significantly  in  the  last 
decade  by  the  rapid  development  of  small  computers,  peripherals  convenient 
for  the  non-specialist,  increased  equipment  reliability,  and  falling  costs. 

A  large  segment  of  the  computer  industry  is  devoted  to  rapidly  improving 
this  area. 

More  specifically,  for  data  analyses  of  the  kind  discussed  in  this 
report,  a  suitable  system  might  consist  of  a  mioroprocessor-based  computer, 
a  floppy  or  hard-disk  driver  system,  a  dot-matrix  impact  printer  (with 
graphics  capability),  a  video  terminal,  and  a  link  to  a  large,  fast  oomr- 
outer.  A  reasonably  satisfactory  system  could  be  obtained  at  a  cost  in 
the  neighborhood  of  ten  thousand  1981  dollars.  This  cost  and  associated 
maintenance  costs  are  small  compared  to  probable  manpower  costs  —  espec¬ 
ially  those  needed  for  the  generation  of  software. 

C.  Computer  Languages.  One  useful  feature  of  the  decentralized,  dedicated 
computer  is  that  it  can  be  made  available  to  one  user  at  a  time.  This 
permits  the  use  of  an  interactive  language  such  as  Forth  (6). 

TJhile  present  versions  of  Forth  are  not  highly  suited  to  involved  num¬ 
erical  analysis,  a  few  remarks  concerning  this  interactive  language  follow 
inasmuch  as  it  seems  likely  that  this  language  and/or  other  interactive 
languages  will  be  developed  in  the  near  future  in  a  way  that  will  greatly 
increase  the  convenience  of  data  analyses  such  as  those  described  herein 
At  present.  Forth  is  highly  suited  to  on-line  use  with  systems  of  instru¬ 
ments  . 

A  Forth-language  resident  system  occupies  about  10K  to  15K  bytes  of 
memory.  It  includes  a  high-level  language,  an  assembly  language,  a  com¬ 
piler,  two  interpreters  (text  and  address),  a  set  of  named  (and  compiled) 
modules,  a  dictionary  containing  the  module  names,  their  addresses,  and 
addresses  which  link  them  to  other  related  modules.  Some  versions  also  have 
text-editing  and  file-handling  software. 
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The  user  "programs"  by  creating  new  named  modules  based  upon  names 
existing  in  the  dictionary,  by  creating  new  named  modules  from  assembly 
language,  by  creating  new  classes  of  named  modules,  including  control 
structures,  data  structures,  and  other  programming  tools.  A  module  is 
compiled  at  the  time  that  it  is  entered  into  the  dictionary  via  the 
keyboard  and  the  new  name  oan  be  used  in  the  same  manner  as  the  original  names 
Thus,  this  language  has  enormous  flexibility  and  extensibility;  the  user, 
in  effect,  creates  a  language  suited  to  the  problem  at  hand. 

The  names  in  the  dictionary  constitute  the  high-level  language.  The 
names  (operators)  are  executed  singly  or  in  sequences  immediately  after 
entry  (in  inverse  Polish  Format)  to  produce  computations  or  input/ output 
(I/O).  Thus,  if  one  withes  to  produce  a  long  sequence  of  computations 
leading  to  a  final  set  of  numbers  and/or  graphs,  an  operator  or  operators 
may  be  exeouted  to  produce  a  segment  of  the  total  with  intermediate  diag¬ 
nostic  I/O  —  or  the  entire  sequence. 

Computer  languages  of  this  kind  are  sometimes  referred  to  as  block- 
diagram  languages  with  the  compiled  modules  being  the  blocks.  The  interp¬ 
reters  parse  the  entered  statements,  consult  the  dictionary,  and  "point" 
the  computer  to  appropriate  modules.  Programs  utilizing  this  kind  of 
software  are  sometimes  referred  to  as  threaded  code,  with  the  addresses 
in  the  dictionary  being  the  thread  that  links  the  modules  at  execution  time. 

The  use  of  a  language  of  this  kind  obviates  the  necessity  for  re¬ 
compiling  and  linking  programs  when  a  small  change  is  made.  Changes  can 
be  made  by  redefining  operators  rather  than  by  changing  programs  —  if  the 
user  prefers. 

While  computers  are  neoessary  for  complex  problems,  they  tend  to 
enslave  the  people  who  interact  with  them.  There  is  a  need  to  devise 
better  modes  of  interaction  to  reduce  the  tedium,  enslavement,  and  dep¬ 
endence  on  large-system  protocol.  Hopefully,  the  foregoing  remarks  can 
be  of  help  in  this  area  of  endeavor.  We  may  look  forward  to  more  power¬ 
ful  computers  and  peripheral  units  at  prices  that  permit  interactive 
computing.  It  also  seems  likely  that  eventually  large  libraries  of  modules 
for  interactive  languages  will  be  developed,  stored  on  floppy  disks,  trans¬ 
mitted  from  one  computer  memory  to  another,  and  distributed  commercially 
at  reasonable  prices. 
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APPENDIX 


N3L  ASC7  COMPUTER  CATALOG  STATUS  TABLE  FOR  NODE  USERCAT/D66/BP0/DUNIIX1/NHP 


SONS:  (  aa 

of  1  August, 

1980) 

HLP1 

SL1DATA5 

0L8YMULT 

HLP3 

0L1DATA5 

LM8YMULT 

HLP4 

LM1DATA5 

SL9YMAT 

HLP5 

IX2HH0S5 

0L9YMAT 

HLP7 

0X1HH0S5 

LM9YMAT 

HLP8 

SL8TUULT 

UJHH0S5 

IN2TGT7  MA2TGT7 


